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Abstract. We consider the focusing nonlinear Schrodinger equations idtu + 
Au + it|u| p_1 = in dimension 1 < TV < 5 and for slightly L 2 supercritical 
nonlinearities p c < p < (1 + s)p c with p c = 1 + -A- and < e -C 1. We prove the 
ON I existence and stability in the energy space H 1 of a self similar finite time blow 

■ up dynamics and provide a qualitative description of the singularity formation 

near the blow up time 



CN . 1. Introduction 
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1.1. Setting of the problem. We consider in this paper the nonlinear Schrodinger 
equation 

iut = - An - V (t, x) £ [0, T) x 



< 



in dimension 1 < N < 5 with 



7V + 2 

1 < p < +oo for N = 1.2 and 1 < p < for N > 3. 

N — 2 

From a result of Ginibre and Velo [7], (jl.ip is locally well-posed in H 1 = H l (R N ) 
and thus, for uq G iJ , there exists < T < +oo and a unique solution u(i) G 
C([0, T), i/ 1 ) to (jl.ip and either T = +oo, we say the solution is global, or T < +oo 
and then lim^T \ Vu(t)\ L 2 = +oo, we say the solution blows up in finite time. 
(jl.ip admits the following conservation laws in the energy space H : 

Q\\ L 2 — norm : f \u(t, x)\ 2 dx = J \uo(x)\ 2 dx; 

Energy: E(u(t,x)) = \ f \ Vu(t, x)\ 2 dx - / \u(t, x)\ p+1 dx = E(u ). 



Momentum : Im(fVu(t,x)u(t,x)dx) = Im{ J X7uo(x)uo(x)dx) 



The scaling symmetry \p- 1 u(X 2 t, Xx) leaves the homogeneous Sobolev space H (Jc 
invariant with 

N 2 , N 

" = 2-J=l- (L2) 
From the conservation of the energy and the L 2 norm, the equation is subcritical 
for a c < and all H 1 solutions are global and bounded in H . The smallest power 
for which blow up may occur is 

4 

which corresponds to a c = and is referred to as the L 2 critical case. The case 
< a c < 1 is the 1? super critical and H 1 subcritical case. 

Even though the existence of finite time blow up dynamics for a c > has been 
known since the 60' using standard global obstructive arguments based on the virial 

l 
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identity, [26J, the explicit description of the singularity formation and of the different 
possible regimes is mostly open. 

1.2. The L 2 critical case and the stable log-log blow up. In the L 2 critical case 
p = p c , a specific blow up regime of "log-log" type has been exhibited by Perelman 
[15] in dimension N = 1 and further extensively studied by Merle and Raphael in 
the series of papers [12], Q3], [32] , Q3], [H], Q2I] where a complete description of 
this stable blow up dynamics is given together with sharp classification results in 
dimension N < 5. The ground solitary wave solution to (II. lj) which is the unique 
nonnegative radially symmetric solution to 

AQ p -Q p + QP = 0, Q p £ if 1 , (1.3) 

see [BJ, [TT], plays a distinguished role in the analysis as it provides the blow up 
profile of log-log solutions. 

Let A be the generator of the scaling symmetry given by A/ = + yVf and 

recall the following Spectral Property which was proved in |12j in dimension N = 1 
and checked numerically in dimensions N < 5 in [3]: 

Spectral Property Let N < 5 and p = p c . Consider the two real Schrddinger 
operators 

A = -A + |pQi + lW~VvQ Pc , C 2 = -A + ^Qr l yVQ Pc , (1.4) 

and the real valued quadratic form for e = e% + i£ 2 G 1 •' 

H(e,e) = (£i£i,ei) + (£2^2, e 2 ). (1.5) 

T/ien i/iere easisis a universal constant 8\ > smc/i £/iai Ve G ii , «/ (ei,Q Pe ) = 
(ei, AQpJ = (e\,yQ Pc ) = (e 2 ,AQ Pc ) = (e 2 ,A 2 Q Pc ) = (e 2 ,VQ Pc ) = 0, then 



H{e,e)>8 1 (J |Ve| 2 + J 



2 P -Ivl 



e e 



We then have the following: 

Theorem 1.1 (Existence of a stable log-log regime, [32], [33], [33], [IS], [E], [IB] . 
[4]). Let N < 5 and p = p c . There exists a universal constant a* > s?<c/i i/iai i/ie 
following holds true. For any initial data uq G if 1 small super- critical mass 

\Q P \ L 2 < \u \ L 2 < \Q P \ L 2 + a* (1.6) 

anc? nonpositive Hamiltonian E{uq) < 0, i/ie corresponding solution to $1.1]) blows 
up in finite time < T < +00 according to the following blowup dynamics: there 
exist geometrical parameters (A(i), x(t), 7(i)) £ R+ x 1^ x R and an asymptotic 
residual profile u* G L 2 st/c/i f/jaf; 

«(*) - -lt—Qp ( m^)*™ - in L? . 

The blowup point converges at blowup time: 

x(t) a?(T) £ R^ as t — >• T, 
i/ie blowup speed is given by the log-log law 



m / iog|iog(r-t)| ^ ^ as t ^ T; (L7) 
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and the residual profile satisfies: 

u* e L 2 but u* i L p , Vp > 2. 

More generally, the set of initial data satisfying hi. 6^ and such that the corresponding 
solution to (DUP blows up in finite time with the log-log law hi. 7]) is open in H 1 . 

In other words, the stable log-log regime corresponds to an almost self similar 
regime where the blow up solution splits into a singular part with a universal blow 
up speed and a universal blow up profile given by the exact ground state Q Pc , and 
a regular part which remains in the critical space u* £ I? but looses any regularity 
above scaling u* £ H a for a > 0. 

1.3. On the super critical problem. The explicit description of blow up dynam- 
ics in the super critical setting is mostly open. In fact, the only rigorous description 
of a blow up dynamics in a super critical setting is for p = 5 in any dimension 
iV > 2, see Raphael [20j for N = 2, Raphael, Szeftel [2T] for N > 3. Note that 
this includes energy super critical problems. In this setting, the existence and radial 
stability of self similar solutions blowing up on an asymptotic blow up sphere -and 
not a blow up point- is proved. These solutions reproduce on the blow up sphere 
the one dimensional quintic and hence L 2 critical blow up dynamic and the blow 
up speed is indeed given by the log log law (II. 7|) . 

In the formal and numerical work [3j, Fibich, Gavish and Wang propose in dimen- 
sion iV = 2, 3 and for p c < p < 5 a generalization of the standing ring blow up 
solutions and investigate a blow up dynamic where the solution concentrates on 
spheres with radius collapsing to zero. Such dynamics are clearly exhibited numer- 
ically and seem to be stable by radial perturbations. A striking feature moreover is 
that these solutions form a Dirac mass in L 2 like in the L 2 critical case: 

M 2 M5 X=0 + \u*\ 2 (1.8) 

for some universal quantum of mass M > 0, and with a specific predicted blow up 
speed: 

\Vu(t)\ L2 a= 5 ~ P (1.9) 

(T-t)i+= (JV-l)(p-l) 

The rigorous derivation of such collapsing ring solutions is an important open prob- 
lem in the field. 



While ring solutions display a stability with respect to radial perturbations, they 
are widely believed to be unstable by non radial perturbations, [3]. In fact, it has 
long been conjectured according to numerical simulations, see [25 an d references 
therein, that the generic blow up dynamics in the super critical setting -at least for 
p near p c - should be of self similar type: 

«(*,*) L^pffEZ^W) 

for some blow up point xt £ and where A is in the self similar regime: 

\(t) ~ Vr^t. 

The delicate issue here is the profile P which does not seem to be given by the 
ground state solution Q p to (II. 3h anymore. In fact, explicit self similar solutions 
may be computed using the following standard procedure. Let explicitly 

1 ( X — Xt\ log(T-t) , 

u(t,x ) = -Q b l^—L\ e —, b>0, X(t) = y/2b(T-t), 

[\(t)]p=i V Mi) ) 
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then u solves (jl.lj) if and only if Qb tP = Qb satisfies the nonlinear stationary elliptic 
PDE: 

AQ b - Q b + ibAQ b + QblQ^- 1 = 0. (1.10) 

Exact zero energy solutions to (|1.10j) have been exhibited by Koppel and Landman, 
[9], for slightly super critical exponent p c < p < (1 + e)p c , eCl, using geometrical 
ODE techniques. The existence of such solutions is related to a nonlinear eigenvalue 
problem and for all p £ (p c , (1 + e)p c ), a unique value b(p) is found such that (jl.lOp 
admits a zero energy solution. An asymptotic law is derived which confirms previous 
formal computations, see [24] pl47 and references therein: 

a c = e ~«5) (1+o(1)) as p -> p c . (1.11) 

Moreover, the self similar profile converges to Q Pc locally: 

Qbip) ~> Q P c in H loc as P Pc- 

However the self similar solutions belong to H 1 n L p+1 but always miss the critical 
Sobolev space due to a logarithmic divergence at infinity: 

\Qb( P )\(y) r w — and hence Qb{p) £ H ac , 

and thus they in particular miss I? and hence the physically relevant space H . 
Eventually, the construction of the self similar solution is delicate enough that it is 
not clear at all how this object should generate a stable self similar blow up dynamics 
for the time dependent problem 



1.4. Statement of the result. The law for the nonlinear eigenvalue (jl.ll)) is 
deeply related to the log-log law (11. 7j) which can be rewritten in the following form: 

| = ^, b = -^, b s = -e-T^)) as S ^+oo. 

This intimate connection between the log-log law and the nonlinear eigenvalue prob- 
lem underlying the self similar equation is at the heart of formal heuristics which 
first predicted (jl.lip . we refer to the monograph [24] for a complete introduction to 
the history of the problem. 

Our main claim in this paper is that the log-log analysis p2] which allowed Merle 
and Raphael to derive the sharp log-log law for a large class of initial data provides 
a framework to somehow bifurcate from the critical value p = p c and prove the 
existence of a stable self similar regime in the energy space for slightly super critical 
exponents, with a blow up speed asymptotically satisfying the nonlinear eigenvalue 
relation (jl.ll)) . In particular, our strategy completely avoids the delicate issue of 
the existence of exact self similar solutions to (jl.lOj) . and provides a framework to 
directly prove in a dynamical way the existence of stable self similar blow up dynam- 
ics, which we expect will apply to a large class of problems. In fact, we will show 
that it is enough to construct a crude and compactly supported approximation of 
the self similar profile (jl.lOj) . Then, the intuition and technical tools inherited from 
the log-log analysis will allow us to obtain new rigidity properties and a dynamical 
trapping of the self similar regime. The outcome is a surprisingly robust proof of 
the existence of a H l stable self similar blow up regime. 
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Theorem 1.2 (Existence and stability of a self similar blow up regime). Let 1 < 

N < 5. There exists p* > p c such that for all p G (p c ,P*), there exists 5(p) > with 
5(p) — ► as p — > pc, i/iere exists > 

<t c = e -^)( 1+ W) (1.12) 

and an open set O in H 1 of initial data such that the following holds true. Let 
uq G O , then the corresponding solution to (jl.lj) blows up in finite time < 
T < +oo according to the following dynamics: there exist geometrical parameters 
(X(t),x(t),j(t)) G R+ x R N x R and an excess of mass e(t) G H 1 such that: 

Vi€[0,T), u (t,x) = ^—[Q p + E(t)](^P-)^ (1.13) 

with 

\Ve(t)\ L2 < 5(p). (1.14) 
The blowup point converges at blowup time: 

x(t) x(T) G R^ as t —>T, (1.15) 
and i/ie Wow t<p speed is self similar: 



vt€[o,r), (i -<%))V2^(p7<^L <(! + %)) v^Cp)- 



(1.16) 



Moreover, there holds the strong convergence: 

Vo-G[0,a c ), u(t)^u* tn as t — > T, (1.17) 

and i/ie asymptotic profile u* displays a singular behavior on the blow up point: 
3R(uq),C(p) > such that: 

Vi?G(0,i?(u o )), C{p)(l-S{p))<-^- ! \u*\ 2 <C(p)(l + 5( P )). (1.18) 

In particular: 

u* G H a for < a < a c but u* $ H aa . 
Comments on Theorem ] 

1. Asymptotic dynamics: Let us stress onto the fact that (|1.13j) does not give 
a sharp asymptotic on the singularity formation, and in particular the question of 
the possible asymptotic stability in renormalized variables in open. In fact, we will 
construct a rough approximate self similar profile b — > Qt, and show that the solution 
decomposes into 

. . 1 r , / x — x(t) y 

u{t,x) = — 2 [Q b(t) +e\ 

with — A(A ~ b(t) and |e(t)|^i <C 1. The key will be to prove a dynamical trapping 
on the projection parameter b(t) and a uniform control on the radiation e: 

Vi€[0,r), (l-5(p))b*(p)<b(t)<(l + 8(p))b*(p) and \e(t)\ 6l « 1, 

but this does not exclude possible oscillations of both b(t) and e(t). In some sense, 
this confirms the analysis in [12], [13] . see also Rodnianski, Sterbenz [23], Raphael, 
Rodnianski [22], where the key observation was that is is not necessary to obtain a 
complete description of the dispersive structure of the problem in renormalized vari- 
ables to prove finite time blow up. The main difference however with these works is 
that we are here in a situation where there holds no a priori orbital stability bound 
on s(t), and a spectacular feature is that rough profiles are enough to capture the 



e »7(*) 
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self similar blow up speed. This somehow confirms the robustness of the log-log 
analysis developed in [13], [H], [15]. Our perturbative approach and the strategy 
of bifurcating from the critical case is also reminiscent from the general framework 
developed by Fibich and Papanicolaou, [5]. 

2. On the behavior of the critical norm: In [T7], Merle and Raphael showed in 
the range of parameters N > 3, < o c < 1, that any radially symmetric finite time 
blow up solution in H 1 must leave the critical space at blow up time with a lower 
bound: 

\u(t)\ A „ e > |log(T - t)\ a ^ as t - T. (1.19) 

The self similar solutions constructed from Theorem 11.21 satisfy a logarithmic upper 
bound -see Remark 14. lb : 

K*)ljf«rc < |log(T-t)|i^T as t^T (1.20) 

which proves the sharpness in the logarithmic scale of the lower bound (I1.19p . A 
logarithmic lower bound could also be derived after some extra work. The singu- 
larity (I1.18P in fact sharpens for this specific class of initial data the divergence 
(|1.20l) . and is according to (j!.18|) the major difference between the L 2 critical blow 
up where the 1? conservation law forces the radiation to remain in the critical L 2 
space, and the super critical blow up where the radiation leaves asymptotically H° c . 

3. Weak blow up: On the contrary to the ring solutions which concentrate in 
L 2 at blow up time according to (|1.8p . the self similar solutions constructed from 
Theorem 11.21 do not concentrate mass according to (jl.l7|) . a situation which is 
referred to as weak blow up, see [21]. A striking feature of the analysis is that the 
L 2 conservation law, even though below scaling, plays a fundamental role in the 
proof of the stabilization of the self similar blow up. 

1.5. Strategy of the proof. Let us briefly state the main steps of the proof of 
Theorem [Ol 

step 1 Construction of an approximate self similar profile. 

Let a small parameter b > and recall that the equation for self similar profiles 
(jl.lQj) does not admit solutions in H . The first step is to construct an approximate 
solution Qt, which is essentially compactly supported and satisfies an approximate 
self similar equation: 

AQ 6 - Q b + ibAQ b + Q b \Q b \ p ~ l = ^ b + 0(a 2 ), 

see Proposition 12.61 This profile incorporates the leading order 0(a c ) deformation 
with respect to the L 2 critical profiles constructed in [13], [IS], while the error ^f b 
is a far away localized error which is inherited from the space localization of the 
profile to avoid the slowly decaying tails. 

step 2 Dynamical trapping of b. 

We now chose initial data such that on a short time, the solution admits a de- 
composition 

u(t,x) = ^—(Q m +s)(t,^P-)e i ^ with |Ve(t)|| 2 «e-^) 
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and aim at deriving a dynamical trapping for the parameter b(t) and the deformation 
e(t). More precisely, introducing the global rescaled time ^ = -h, the dynamical 
system driving A is 

-AA t = -y ~b(t) 

and hence finite time blow up in the self similar regime will follow from 

~ 6 > 0, |Ve(t)||a <e~^ (1 ~ c) 

for some small constant c > in the maximum time interval of existence, see Propo- 
sition [27121 In order to derive such dynamical controls, we run the log-log analysis 
developed in [il5j by keeping track of the leading order 0(a c ) deformation. The out- 
come is the derivation of two structural monotonicity formulae for the parameter 
b. The first one is inherited from the local virial control first derived in jT2j and 
roughly leads to 

<J C + |Ve|| 2 -e - T < b s , 
see Proposition 13.31 The positive term +a c in the above LHS is a non trivial 
supercritical effect and is a consequence of the structure of the self similar profiles 
in the local range \y\ < 1. The second monotonicity formula is a consequence of the 
L 2 conservation law and the control of the mass ejection phenomenon: 

b s < cr c - e~^, 

7T 

where the nonpositive term — e" is obtained from a flux computation in the far 
away zone \y\ » 1 which is based on the presence of the slow decaying tails of self 
similar solutions. The outcome is the derivation of the dynamical system for b: 

b s ~ a c - e~b 

which traps b around the value 

b ~ b* with <7 C ~ e~b*. 

Note that this shows that the self similar blow up speed is derived from the con- 
straints both on compact sets and at infinity where the dispersive mass ejection 
process is submitted to the global constraint of the L 2 conservation law. Note also 
that the L 2 critical log-log law corresponds to the dynamical system 

7 -21 

b s ~ — e b 

and hence the supercritical self similar blow up appears as directly branching from 
the L 2 critical a c = degenerate log-log blow up. The nontrivial pointwise control 
on e 

|Ve(t)|£ a « e~ (1 ~ c) ^) 

will also follow from the obtained Lyapounov controls. Here a difficulty will oc- 
cur with respect to the L 2 critical case to control the nonlinear terms in e due in 
particular to the unboundedness of the scaling invariant critical Sobolev norm. A 
new strategy inspired from [2T], [23J, |22J is derived which relies on the control of 
Sobolev norms strictly above scaling in the self similar regime, see section 14.11 
The conclusions of Theorem 2 are now a simple consequence of this dynamical trap- 
ping of the solution. 

This paper is organized as follows. In section [2l we construct approximate self 
similar solutions, Proposition E21 We then describe the set of initial data leading to 
self similar blow up, Definition 12.81 and set up the bootstrap argument, Proposition 
12.121 In section [U we derive the key dynamical controls and the two monotonicity 
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formulae, Proposition 13.31 and Proposition 13.71 In section HJ we close the bootstrap 
argument as a consequence of the obtained Lyapounov type controls and conclude 
the proof of Theorem 11.21 

Acknowledgments. F.M. is supported by ANR Projet Blanc OndeNonLin. P.R 
and J.S are supported by ANR jeunes chercheurs SWAP. 

Notations We let \V\ a be the Fourier multiplier jVp/(£) = \Z\ a f(Q- We let Q p 
be the unique radially symmetric nonnegative solution in H 1 to 

AQ P - Q p + Q p p +1 = 0. 

We introduce the error to 1? criticality: 

4 N 2 N(p- Pc ) 

where < a c < 1 is the Sobolev critical exponent. We introduce the associated 
scaling generators: 

Af = ^—f + yVf, Df = ^f + yVf = Af + a c f. (1.22) 
p—1 2 

We denote the L 2 (W N ) scalar product 

(f,g) = / f{x)g(x)dx 



and observe the integration by parts formula: 

(Df,g) = -(f,Dg), (A/, g) = -(/, Kg + 2a c g). (1.23) 
We let L = (L + ,L-) be the linearized operator close to Q p : 

L + = -A + l- P QP' 1 , L_ = -A + l- QP-\ (1.24) 



2. Description of the blow up set of initial data 

This section is devoted to the description of the open H 1 set O of initial data 
leading to the self similar blow up solutions described by Theorem 11.21 which relies 
on the construction of approximate self similar profiles. 

2.1. Construction of approximate self similar solutions. Our aim in this 
section is to construct suitable approximate solutions to the self similar equation 
(|1.10|) . Let us make the following general ansatz: 



and introduce the rescaled time 

ds _ 1 
dt ~ A*(t)' 

then u is a solution to (jl.lj) if and only if v solves: 

id s v + At; — v — i-^-Av + v \v | p_1 = (7^ — l)v + i— - • Vv. 

A A 



Let us fix 
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and look for solutions of the form: 

v(s,y) = Q b (s)(y) 
where the unknowns are the mappings 

s^b(s), b^Q b . 

This corresponds a to a slow variable formulation of a generalized self similar equa- 
tion which goes back to previous formal works, see [24], and was rigorously used 
in [IB], p3], see also [TU| . [22] for related transformations in different settings. To 
prepare the computation, we let 

b s = o-cVb 

where /J, b is a function of b which will be made explicit later on. The generalized 
self similar equation becomes: 

ia ^^ + A ® b + ibA ® b + QbW' 1 = o. (2.i) 

Let us perform the conformal change of variables 

ib\y\ 2 

Pb = Qbe 4 , 

then a simple algebra leads to: 

OP 1 

ia ^ b ^f + APb ~ Pb ~ i(TchPb + 4^ 2 + a cVb)\y\ 2 Pb + PbW 1 = 0. (2.2) 

Our aim is to find fi b so as to be able to construct an approximate solution to (|2.2p 
with an error of order formally a 2 in the region |y| < j > . Our construction is ele- 
mentary and relies on the computation of the first term in the Taylor expansion of 
(Mb) Qb) m °~c near the L 2 critical value o~ c = 0. 

Let a small parameter < rj « 1 to be fixed later, a non zero number 6, and set 

2 



Rb = jkjVl " V and R£ = y/l- vRb- (2.3) 

Denote B Rb = {y E R N , \y\ < R b } and dB Rb = {y £ R N , \y\ = R b }. We introduce 
a regular radially symmetric cut-off function 

for Id > R b , , , . . . 

1 for N i R^, ° ^ <*(x) < 1, (2.4) 

such that: 

\(f>' b \ L oo + |A0 6 |ioc -»• as |&| — »- 0. (2.5) 

We also consider the norm on radial functions ||/||cj = niaxo<fc<j- ||/^( r )llL°°(R+)- 
Let us start with the construction of the profile for a c = where we view a c and p 
as independent parameters and leave p supercritical. 

Proposition 2.1 (Q b profiles). There exists p* > p c and C,rj* > such that for 
all p c < p < p* , for all < r\ < 7]* , there exists b* {rj) , e* (rj) > going to zero as 
rj — > such that for all \b\ < b*(rj), the following holds true: 

(i) Existence of a unique P b profile: there exists a unique radial solution to 
' Ap (o)_ p (o) + « p (o ) + ^(o)y = 0) 

pf >0 in B Rb , ' (2-6) 

[ P b (0) (0)e(Q p (0) -£*(»?), Q p (0)+e*^ p(°\ Rb ) = 0. 
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Moreover, let 



then P^ is twice differentiable with respect to b 2 with uniform estimate: 



p(o) {r) = P m {r)Mr) 



8(|ft|r) 



P 



(0) 



Qp 



C3 



6(\b\r) I dP, 



(0) 



db 2 



as b 



d 2 (b 2 



< C 



C3 



where 



6(w) 



1 



0<w<2 



61(2), 



— dz + l w>2 



anc? /9 is i/ie unique solution in H^ ad to 



l+p = ^\y\ 2 Q- 



(ii) Properties of the profile: Q^' = e 4 P b w satisfies: 



(0) 



,(0) 



(°)ir)(°)|p-i 



(o) 



,(0)^64- 



(0) - 2V</> 6 • VP fc {0) + P 6 (0) (A^) + - h){P^) p 



and for any polynomial f(y) and integer k = 0,1, 



f(y) 



Op 



for some constant c p depending on p. 

(Hi) Computation of the momentum and the L 2 mass: there holds 



Im 



vgf>g<°> 



0, Im 



(0)|2 
b \2 



and the supercritical mass property: 
d 2 



db 2 



\Q. 



(0)|2 



1 6=0 



(2.7) 

0, (2.8) 
(2.9) 

(2.10) 

(2.11) 

(2.12) 
(2.13) 

(2.14) 



(2.15) 



cq(p) with cq(p) -> co(p c ) > as p -> p c . (2.16) 



The proof of Proposition 12.11 is parallel to the one of Proposition in [13] and 
Proposition [14] and relies on standard elliptic techniques and the knowledge of the 
kernel of the linearized operator close to Q, explicitly: 

Ker(L + ) = span(VQ), Ker(L_) = span(Q), (2.17) 

see [25], [2], and the fact that we are working here before the turning point | and 
hence with uniformly elliptic operators. The detailed proof is left to the reader. 

Remark 2.2. The question of the value of the energy of the modified profile is an 
important issue. It will indeed be computed in Proposition \2. 6\ for the full approxi- 
mate profile, see (|2.36p . 

After the turning point |, leading order phenomenons are of linear type and a 
natural prolongation of the approximate blow up profile is given by the so called 
outgoing radiation. 
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Moreover, let 9 be given by $2.10)) . and consider 

T h = lim \y\ N \C b (y)\ 2 , 



lim \y\ 

\y\-y+oo 



then there holds: 



\y\-(\( b \ + \y\\v(C b )\) 



L°°(\y\>R b ) 



< +00, 



Lemma 2.3 (Linear outgoing radiation). See Lemma 15 in |14j . There exist uni- 
versal constants C > and r/* > such that VO < rj < rf , there exists b*(rj) > 
such that\/0 < b < b*(r]), the following holds true: let be given by 112.12]) . there 
exists a unique radial solution Cb t° 

'.r,>.2 ( 2 - 18 ) 

|vc 6 | 2 < +00. 

(2.19) 

(2.20) 
(2.21) 

(2.22) 

(2.23) 
36** (77) 

(2.24) 
(2.25) 



IvCbl 2 < r 



1-Ct] 



For \y\ large, we have more precisely: 
V|y| > R 



? 2 



e -2(l-C,)ffii > | yr|C6(y)| 2 



> -r 6 > e 

5 



V\y\>R 2 b , \VCb{y)\< 



c 



r? 



iV 
' 2 



|yr 1 2 I 6 ' 

For \y\ small, we have: \/a G (0,5), 3?7**((j) smc/i that VO < rj < rj**(a) 
such that VO < b < b**(rj), there holds: 



Cb(y)e-^ 



<r 6 2 



C 2 (\y\<R b ) 

Last, C& is differentiable with respect to b with estimate 



Ob 



f(y> 



dy k 



Remark 2.4. Recall from (|2.10|) that 6(2) = j. Moreover it is enough for our 
analysis to compute the radiation with the L 2 scaling generator D in (12 . 18|) instead 
of A as the error will generate lower order terms. 

Remark 2.5. Now that Tj, has been defined in (|2.19p . we can give a more precise 

formulation of the estimate (|2.14|) satisfied by which is a direct consequence 
of (|2.8|) and the localization procedure (|2.13p .- for any polynomial f(y) and integer 
fc = 0,l, 

8il '" <ri"- C ". (2.26) 

L°° 

The proof of Lemma [2~73l is completely similar to the one of Lemma 15 in [14] -see 
also Appendix A in [15]- and hence left to the reader. In particular, the nondegen- 

eracy (j2.22|) is a consequence of the nonlinear construction of the profile Pff^ . 

One should think of as being an approximate solution to the generalized self 
similar equation (|2.1|) with /x& = -self similar law- and an error of order a c in the 
elliptic zone |y| < ^. We now claim that there exists a -locally unique- non trivial 
fJ>b(p) > which allows one to construct an approximate solution to the generalized 
self similar equation (12.11) of order O(o~ 2 + V 2 ) in the elliptic zone \y\ < |: 
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Proposition 2.6 (Approximate generalized self similar profiles). There exists p* > 
p c and C,ci,r]* > such that for all p c < p < p* and < 7] < r]* , there exists 
co(p) > 0, 6*(ry) > such that for all \b\ < b*{rj), the following holds true: 
(i) Construction of the modified profile: there exist \x b = fJ>(b,p) > and a radially 
symmetric complex valued function T b = T(b, p) with 



+ 



db 



+ 



c-' 



dfi b 



db 



<C as b^O, (2.27) 



fJ>b 



pit 2 



(l + 2a c )\yQ p \ 2 L2 
with the following properties. Let 



as b — > 



P b = P h {0) +a c T b , Q b = P b e l 



then 



-KT c fJ, b - 



dQb 
db 



AQb + Qb-ibAQb-QblQ, 



m-l 



*[ 0) + * 



(i) 



with ^[ 0) given by (l2~T2ll satisfies (l2~26l) . and for k = 0, 1: 



dy k 

L°°{\y\<RD 

for some constant c\ > depending only on N and: 



< a. 



l + Cl 



dy k 

L°°(\y\>B-) 

(ii) Estimate of the invariants on Q b : there holds 



<Ca, 



Im / VQ b Q l 



0. 



(/ 



Im\ y- VQbQi 



\Q b 



~|l/Q P |l(l + 0(|6|+ff c )) as b 
= J \Q p \ 2 + M(b) + 0(a c ) 



0, 



with 



d 2 



M(0) = and -^M{b)\ b=0 = c (p) -» c (p c ) > as p 



and the degeneracy of the Hamiltonian: 

\E(Q. 

Proof of Proposition 12.61 



<rJ- c " + c<7 C 



(2.28) 



(2.29) 
(2.30) 

(2.31) 

(2.32) 

(2.33) 
(2.34) 

(2.35) 
(2.36) 



The proof relies on a Taylor expansion in a c of formal solutions to (12. The 
choice of fi b is dictated by the presence of a non trivial kernel for the operator L_ 
driving the imaginary part of (|2.ip near Q p as given by (|2.17l) . while L + on the 
real part is invertible in the radial sector. The computation of the invariants and 
in particular the energy degeneracy t|2.36|) follow from Pohozaev identity. 



step 1 Construction of T b . 
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Let Q b = P b e~ ibh ^~ and ^ b given by (l2~29fl , then: ^ b = * 6 e~ i&1 4~ with: 

dP 1 

- #6 = *VcM6-^ + AA - P 6 - m c 6P 6 + -{b 2 + a c fi b )\y\ 2 P b + AW" 1 - (2-37) 



Let be given by (|2.13|) . equivalently: 



(0) 



AP, 



(0) 5(0) 



|2 6(0) , 6(0) i 6(0) m-1 



4 \y\ r b r b \ r b 



(2.38) 



We expand P& = P^ + a c T b . Let 06 be the cut off function given by (12, 4| . We 
compute: 



1 



1 



96 

3(0) 



+ (P fe W + a c T b )\P b W + o-c^r 1 - (P ( 
+ cr c 



(0) 



-(L fe)+ Pe(T fe ) + ^| y | 2 pW 



6 



^(P^f-^e^) - ia c {P^f-Hm{T h 



B(0)\p-1 



where we introduced the linearized operators close to P, : 



dP 

-{L h )-Im{T h ) + 

(0). 



(0) 



bP 



(o) 



A + 1 - -^| y | 2 - p{Pr) p -\ (L-) b = -A + 1 - -fclyf - (P 6 W ) 



(0)\p-l 



We thus aim at finding {fjL b ,T b ) so as to cancel the 0(a c ) in the RHS of (|2.39|) : 



o p 

(L b )+Pe(T 6 ) + ^|y| 2 P fe (0) = 0, -(L b ).Im(T b ) + 



(o) 



bP b {0) = 0. (2.40) 



Recall that in the limit 6^0, (L+)o = £+ is an elliptic invertible operator 
in the radial sector and (LJ)q = L_ is definite positive on (Span(Q)) with 
Ker(L_) = Span(Q). The following lemma is a standard consequence of the Lax- 
Milgram theorem and the perturbative theory of uniformly elliptic Schrodinger op- 
erators, and its proof is left to the reader: 

Lemma 2.7 (Invertibility of (£+)&, (L^) b ). Given r\ > small, there exists b*(rj) > 
such that for all \b\ < b*{rj), the operator (L + ) b is invertible in the radial sector 
with 

\\e {1 ' v) W (L + )^f\\jfl < CJe^-^ — fh?. (2.41) 

Moreover, (£-)& admits a lowest eigenvalue X b with eigenvector which are C 1 
functions of b with: 



and 



Moreover, 



(!-,)« 



0b 



+ 



b| (Cb-Qp) 

(l_ r;) £(l^)^6 

db 



as 6^0 



B(\b\r) 



<C as 6^0. 



(2.42) 
(2.43) 
(2.44) 



V/ E (Spanfo)) 1 , ||e l 6 l (L^.)^ 1 f\\ H 2 < C n \\e { ^"'^T f\\ L ,. 

From (I2.44D . the solvability of the second equation in (|2.40l) imposes the choice 
of Hb~ 



M6 



06 



&Pi 0) ,6 =0 i.e. M 6 



(0) 



6 



(2.39) 



14 F. MERLE, P. RAPHAEL, AND J. SZEFTEL 

We now observe the crucial non degeneracy of /i b as b — > from (|2.8j) , (|2.42j) and 
(J23H): 

where we used the computation from (|2.1ip and L + (AQ P ) = —2Q P : 



The non degeneracy of the denominator as b — > and the uniform differentiability 
properties (|2.9|) . (|2.43|) ensure that /j, b is a C 1 function of b with 

dfib <C as b^O. (2.46) 



6*6 

Hence from (l2~4T1) . (jOID and the uniform bounds JZS]), (l2~45l) . (fQ6l) we may find 
T b solution to (|2.40|) which is a C 1 function of b satisfying the uniform bounds (|2.27j) , 

step 2 Estimate on the error. 

We now turn to the proof of the estimate of the error (I2.30P which amounts esti- 
mating the remaining terms in the RHS of (|2.39j) . The nonlinear term is estimated 
thanks to the homogeneity estimate: 

(P b (0) + a c T b )\Pl 0) + a.nr 1 - (Pl 0) y - a c p{P^y- l Re{T h ) - ia^P^f^Imin 

a p c \T p \P for p < 2, 
a p c \T p \P + a 2 c \Pl 0) \P- 2 \T p \ 2 for p > 2, 

and PJZZD now yields (HOOt (123111 . 
step 3 Computation of the invariants. 

(|2T32l) still holds because Q b is radially symmetric. (12T331) . (I2~34ll . (I2~35ll follow 

from (|2.15p . (|2.16p . the decomposition Q b = + a c T b e~ lb ^~4~ , and the fact that 
Im(Tfc) = 0(6) by (|2.40|) . To compute the energy, we use the Pohozaev multiplier 
AQ b on (I2.29p . We first integrate by parts to get the general formula: 



< 



Re(AQ b -Q b + ibAQ b + Q b \Q b \ p -\ AQ b ) = -2E(Q b ) + a c (^2E(Q b ) + J \Q 
and hence from (|2.29|) : 



2 

b\ 



2E(Q b ) - Re ( AQ b , tf 6 + ia c p, b dQh 



db 

2E(Q b ) - (Re(* b ), AS) - (Im(* b ), AG) + a cfM) 



d& . \ /as 

— , AS - — , A9 
96 ' J {Ob' 



= a c (2E(Q b ) + J \Q b \ 2 ^j , (2.47) 

where S, are defined by Q b = S + iO. This together with (12.34(1 and the estimates 
on ty b (12361) (I2T301) ([2~3Tjl yields the degeneracy (l2~36]l . 
This concludes the proof of Proposition 12.61 
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2.2. Setting of the bootstrap. We are now on position to describe the set of ini- 
tial data O leading to the self similar blow up and to set up the bootstrap argument. 
Pick a Sobolev exponent 

^{^{iwWj) (2 - 48) 

independent of p and close enough to 0. 

Definition 2.8 (Geometrical description of the set O). Pick a number uq > small 
enough. Then for p £ (Pc,P*{ v o)) with p*{vo) close enough to p c , we let O be the 
set of initial data uq £ H 1 of form: 



\ p-i 
A 

for some (Ao, &o, ^0, 7o) £ R+ x x R N x R with the following controls: 

(i) bo is in the self similar asymptotics fll.llj) : 

C' < a c < Tl o "° ; (2.49) 

(ii) Smallness of the scaling parameter: 

< Ac < If ; (2.50) 
(m,) Degeneracy of the energy and the momentum: 

A o l-^o I + A o 



< lt°; (2.51) 



Im y^J Vuouq 

(iv) H 1 n H a smallness of the excess of L 2 mass: 

j ||Vre | 2 + / |Ve | 2 + / \eo\ 2 e~\y\ <T\^. (2.52) 

Remark 2.9. Hidden in Definition \2.8l is the choice of the parameter 7] > entering 
in the construction of in Proposition \2.1l In all what follows, we will need the 
fact that given a universal constant C > 0, we have the control 

T l-Cr, < pl-^o 50 
1 b — 1 b 

This holds provided r\ > has been chosen small enough with respect to uq. We may 
for example take 

V = -o 100 - 

Now, pushing ij —* requires pushing b — > or equivalently p —* p c from (|2.49j) . 
This is how the asymptotics (|1.12j) follows. 

Remark 2.10. Observe that the set O is a non empty open set in H 1 . Indeed, pick 
a small parameter vq and p close enough to p c . Pick bo > such that (12.491) holds, 
and pick then Ao > such that (12.501) holds. Let f(y) be smooth, real, radial and 
compactly supported in the ball \y\ < 1 and such that (f,Q) = 1, and let Eq = Hof 
with no to be chosen. From 

d -E(Q bo +Hof)\v =o = -(f,Q)(l + o(l)) = -l + o(l) as b - 0, 



d^o 

and the degeneracy of the Hamiltonian (|2.36|) . we may find fio = 0{T l ~ 2vo ) such 
that \E(Q bo +e )| < 1^°° so that (I23T1) /iotas. (12321) now follows from the size of 
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(2.53) 



Let now uo G O and u(t) be the corresponding solution to (II. lh with maximum 
life time interval [0, T), < T < +oo. Using the regularity u G C([0,T), H 1 ) and 
standard modulation theory, we can find a small interval [0,T*) such that for all 
t G [0, T*), u(t) admits a unique geometrical decomposition 

u{t,x) = -j—(Q m+ e) U^@) e**> 
A^(i) V A(t) ) 

where uniqueness follows from the freezing of orthogonality conditions: Vi G [0, T*], 

(ei(t),|y| 2 S) + ( £2 (t),|y| 2 e) =0, (2.54) 

(£i(*),yE) + (e 2 (t),ye) = ) (2.55) 
(e 2 (t),A 2 S) - ( £l (i),A 2 e) =0, (2.56) 
(e 2 (t),AIl)-(ei(t),Ae)=0, (2.57) 

where we have denoted: 

e = ei+ie 2 , Qb = S + iO 

in terms of real and imaginary parts. See [12], [13] for related statements. Moreover, 
the parameters (A(i), b(t),x(t),^(t)) EK;xi;xK w xl are C 1 functions of time 
and e G C([0, T*),!! 1 ) with a priori bounds: Vt G [0,T*), 



r 6 J)°<^<r— °, 
o<A(t)<r 6 1 ° ) , 



Im VuqMo^ 



[A(t)] 2 (i-^)|£b| + [A(t)] 1 - 2 ' 7 ' 

| |Ve(i)| 2 + | k(i)| 2 e-^<r 
/ ||V| CT e(t)| 2 < r 



<r r 10 

<i 6(*)' 



l-20f 

6(t) ' 



1-501/q 
6(t) • 



(2.58) 
(2.59) 

(2.60) 
(2.61) 
(2.62) 



Remark 2.11. T/ie sihci if 1 subcriticality of the problem implies: 

N 2 N N 
~2 



Or 



N 

p-1 < 2 p+1' 



Hence from Sobolev embedding, the H a control (|2.62|) together with (|2.61|) ensures 
for a close enough to a c : 



IP+1 < 



< 



pi— 50^o 
1 b 



2 



< pl+^0 



/or some universal constant zq > independent of p for p close enough to p c 

Our main claim is that the above regime is a trapped regime: 
Proposition 2.12 (Bootstrap). There holds: Vt G [0,T*) ; 

1 b(i) - ac - 1 6(t) ' 



(2.63) 



< A(t) < r 20 



&(*)> 



r 20 



| |Ve(t)| 2 + | | £ (t)| 2 e-^<r 



l-10fo 
6(t) ' 



(2.64) 
(2.65) 

(2.66) 
(2.67) 
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V\ a e{t)\ 2 < F 1 -^ (2.68) 



b(t) 
and hence 

The next section is devoted to the derivation of the key dynamical controls at 
the heart of the proof of the bootstrap Proposition 12.121 which is proved in section 
IH Theorem 11.21 will be a simple consequence of Proposition 12.121 



3. Control of the self similar dynamics 

In this section, we exhibit the two Lyapounov type functionals which will allow 
us to lock the self similar dynamics and prove the bootstrap Proposition 12.121 The 
key is the dynamical lock (|2.64l) of the geometrical parameter b(t) which controls 
the selfsimilarity of blow up from the modulation equation 

b ~ — ^ = — XXt- 
A 

This corresponds to an upper bound and a lower bound of b. The proof will fol- 
low by somehow bifurcating from the log-log analysis and by tracking the effect 
of the leading order a c deformation in the Qb profile. The lower bound on b is 
a consequence of the local virial estimate type of control first derived in jT3j, see 
Proposition 13.31 the upper bound follows from the sharp log-log analysis derived in 
[T5] and uses very strongly the I? conservation law, see Proposition 13.71 Remember 
also that we have no a priori orbital stability bound on neither b or e, and indeed 
the upper pointwise control (|2.67l) on e requires both monotonicity properties, see 
step 1 of the proof of Proposition 12.121 in section 14.11 

3.1. Preliminary estimates on the decomposition. Let us recall the geomet- 
rical decomposition fj2 .53jl : 

u(t,x) = -J—(Q b{t)+ e) U^P) e**«> 

and derive the modulation equations on the geometrical parameters and preliminary 
estimates inherited from the conservation laws. We let the rescaled time 

f* dr 

s (t)= / s* = s(T*) G (0, +oo], 



/o A2(r)' 

and compute the equation of e in terms of real and imaginary parts on [0, s*): 

(b s - Wb)-^ + d s e 1 - M_(e) + bA £l = ( y + b\ AS + 7,6 + y • VS 

+ + bj Aei + %e 2 + y • Vei 

+ Im(^ b )-R 2 (e) (3.1) 

BP) /A \ r 

(b s -a cl i h )— + d s e 2 + M + (e) + bkE 2 = f -f + b) AG - %H + ^ • V6 

+ + bj Ae 2 - %ei + y • Ve 2 

- fle(* 6 )+i2i(e), (3.2) 
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with j(s) = -s+j(s) and * b given by ()2,29jl . Here we also denoted M = (M+, M_) 
the linear operator close to Q^, explicitly: 

M+(e) = -A £l + ei -(l + (p- 1) I W^i - (P - l)se|Q b | p - 3 e 2 , 

M_(e) = -Ae 2 + e 2 - + (p - 1) jji) W" 1 ^ - (p - l)Ee|Q 6 r 3 £l . 
The non linear interaction terms are explicitly: 

fla(e) = (£i + E)| £ + Q 6 r 1 -S|Q b r 1 (3.3) 

- (l + (p - !) ^2) l^r^i - (P - l)ZO\Qb\ p ~ 3 e 2 , 

R 2 (e) = ( £2 + e)| £ + Q fe r 1 -e|Q fe r 1 (3.4) 

- (i + (p- !) ^2) \Qb\ p - x e2 - (p - i)se|Q 6 r 3 £l . 

We now claim the following preliminary estimates on the decomposition: 

Lemma 3.1. There holds for some universal constants C > 0, 5(p) —* as p ^ p c 
and for all s G [0, s*); 

(i) Estimates induced by the conservation of energy and momentum: 

|2(ei,E) + 2(e 2 ,e)| <c(J \Ve\ 2 + J \s\ 2 e~\A + T\- llv ° (3.5) 

|(e 2 ,VS)| <S(p) (J |Ve| 2 + J | £ | 2 e -^V + r^ 50 " . (3.6) 
fii,) Estimates on the modulation parameters: 



+ \b s \<c(J \Ve\ 2 + j | e | 2 e -i^+r 



l-lll/O 

6 ' 



(3.7) 



+ 



1 



< 5(p) (| iVelV 2 ^)^ + /'| £ |VI^' 



+ Cy |Ve| 2 + r£- llMD . (3.8) 

Proof of Lemma 13.11 

It relies as in [15] on the expansion of the conservation laws, the choice of orthog- 
onality conditions for e and the bootstrapped controls (I2.58L (|2.59|) . (|2.60p . (|2.61|) . 
(I2~62l) . 

stepl Expansion of the conservation laws. 



(|3.5|) and (13. 6p follow from the expansion of the momentum and the energy using 
the decomposition (|2.53p and the estimates of Proposition 12.61 
For the momentum: 

2Im(e,VQ b ) = Im J Vee - A 1_2o " c Im (J Vu wj 
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from which using fl2 .60|) : 

\(e 2 ,VQ)\ < 5{p)(^j\Ve\ 2 + j\e\ 2 e-\y\y +C\e\ 2 tih +Tf 

< S{p) (J \Ve\ 2 + f \e\\-\v\y +T 1 -^ 

where we interpolated between H a and H 1 and used (|2.6ip . (|2.62p . 
For the energy: 

2 ( Sl , S + bAQ - Re(V b ) + Wh^) + 2 (e 2 , 6 - 6AS - Im(%) - 
= -2X^ 1 -^E + 2E(Q b ) + (M + (e),e 1 ) + (M^(e) 1 e 2 )- J \e\ 2 

F(e) (3.9) 



p + 1 

where -F(e) is the formally cubic part of the potential energy: 



F(e) = \Q b + e\ p+l - \Q b \ p+1 - (p + l)ite(e, Q^p 1 ) 

- (p + 1) (l + (p - 1)^) IQal^e? - (p + 1) (l + (p - |Q 

- 2(p + l)(p-l)m\Q b \P~ 3 e 2 e 1 . (3.10) 



6| E 2 



We then use standard homogeneous estimates and Sobolev embeddings like in [12] . 
[T3] to estimate the nonlinear term F(e). Indeed, choose \i such that < fj, < 
min(l,p — 1), then by homogeneity: 

\F(e)\ <\s\p +1 + \Q b r l ^\e\ 2+ ^ 

and thus from Holder with ^7 = § + ^py and the bootstrap controls (|2.61l) . (|2.63|) : 



1^001 < \e\ p+l 



a(2+fM) , ,(l-a)(2+M) 



L 2 



< r^ + frj- 20 ^^^^^ (3.11 



for some zq > independent of p. Injecting this into (13.9(1 together with the de- 
generacy estimate (|2.36l) . the bootstrap bound (j2.58|) . (I2.59P and the orthogonality 
condition ([2~57]) yields (p£S)) . 



Remark 3.2. Note that the algebraic formula (13. 9p e?oes noi use the orthogonality 
conditions on e. 

step 2 Computation of the modulation parameters. 

The estimates (|3.7I) . (|3.8I) are obtained by computing the geometrical parameters 
from the choice of orthogonality conditions (|2.54l) . (|2.55j) . (|2.56l) . (j2.57j) and by 
relying on the estimates induced by the conservation of energy and momentum 
(|3.5p (13.6(1 . The computation is the same like the one performed in [15] up to 
0(a c ) = 0{p — p c ) terms which are brutally estimated in absolute value using the 
bootstrap bound (j2.58|) . The detail of this is left to the reader. 
This concludes the proof of Lemma IBTTl 
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3.2. Local virial identity. We now proceed through the derivation of the local 
virial control. The corresponding monotonicity property was first discovered in |12j 
and will yield a strict lower bound on b. 

Proposition 3.3 (Local virial identity). There holds for some universal constant 
c\ > the lower bound: 

Vs 6 [0, s*), b s > ci (a c + J | Ve| 2 + J \e\ 2 e~^ - T 1 '^ . (3.12) 

Remark 3.4. The super critical effect relies in the presence of the term a c in the 
RHS of (|3.12|) . The positive sign is crucial and structural and is the reason why we 
had to push the construction of Q b to an order a 2 . 

Proof of Proposition 13.31 

stepl Algebraic derivation of the b s law. 

The first step is a careful derivation of the modulation equations for b. For fur- 
ther use, we shall exhibit first a general formula which does not rely on the specific 
choice of orthogonality conditions in e. Once the formula is derived, the use of the 
orthogonality conditions and suitable coercivity properties inherited from the Spec- 
tral Property stated in the introduction will yield the claim. 

Take the inner product of (|3.ip with (— AG) and (|3.2I) with AE and sum the 
obtained identities using (|1.23h to get: 

((f - ae ) - (f ' Ae ) - (-> A f ) + + - - 

- (M + (e) + 6Ae 2 ,AS)-(M_(e)-6Ae 1 ,Ae)-7 s {(ei,AS) + (e 2 ,AG)} (3.13) 
y + bj {(e 2 , A 2 £ + 2<r c A£) - ( Bl , A 2 G + 2<r c AG)} - y • {(e 2 , VAS) - ( £l , VAQ)} 

- Re (^AQ b , * 6 + ia^^j + a c %\Q b \ 2 L , + {R^e), AS) + (R 2 (e), AG). 

We compute the linear term in e in (|3.13l) . For this, consider the Q b equation (|2.29l) . 

2 

compute the equation satisfied by Q(fiy) and differentiate the obtained identity 
with respect to fj, at (J, = 1. This yields: 



M + (AQ b ) + 6A 2 Q 



M_(AQft) - 6A 2 S 



£ + 6AG - Re(%) + <r c y. b -^ 



G - 6AS - Im(V b ) - OcHb-r^ 



dG 

+ Re(A^ b )-a c fi b A—. 



+ Im(A^ b ) +a cf i b A—. 
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Integrating by parts and using the conservation of the energy (|3.9I) . we obtain the 
following identity: 

-(M+(e) +6Ae 2 ,AS) - (Af_(e) - 6Ae x ,AG) 
= -(ei, M+(AQ 6 ) + 6A 2 Q + 26ct c AG) - (e 2 , M_(AQ 6 ) - frA 2 S - 26ct c S) 



£i,2 



e 2 ,2 



£ + 6AG - Re(^ b ) + cr c /i b 



99" 
~db 
<9S 

G - 6AS - Im(¥&) - ^cMftflT 
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i?e(A* 6 ) + a c /i 6 A— - 26<r c AG 
9S 

Im(A¥&) - a c ii h k— + 2fcr c AS 



2fcr c [-(ei, AG) + (e a , AS)] - (e 1} lte(Atf 6 )) - (e 2 , im(A¥ 6 )) 



db 



ei,A— - £ 2 ,A— 



as 
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+ (M + ( £ ), £l ) + (M_( £ ),£ 2 )- / |e| 2 

J P + 1 

= <y(p)<r c - (si, ite(A* 6 )) - (£ 2 , Im(A* 6 )) - 2A 2 ( 1 -^£ + 2£(Q 6 ) 
+ (M+(e),ei) + (M_(e),e 2 )- /|e' 2 



2A 2 ( 1 -^) J B + 2£(Q 6 ) 
2 ^(£) 



p+ 1 



F( £ ) 



where we recall that <5(p) denotes a generic constant 5{p) — > as p — ► p c . We now 
inject this into (|3.13|) , A key here is to use the Pohozaev identity ((2.4TB : 



2E(Q b ) - Re ^AQ b , V b + i^b^f) = o c \2E{Q h ) + J \Q b \ 2 ^j 



to generate a nonnegative term in the RHS of (13 . 13B : 



de 



dS 



3S\ 



bs -or, AS - — , AG - e 2 , A— + Bl> A— + {(e a , AS) - (ei, AG)} 



db J 



db J 



db' J \db 
5(p)a c + a c {2E{Q b ) + \Q b \\ 2 + %\Q b \ 2 L2 ) - 2X 2 ^-^E 

- (£ 1 ,i?e(A* 6 ))-( £2 ,Jm(A* 6 )) 

+ (M + (e),e 1 ) + (M^e),e 2 )- J \e\ 2 + (R^e), AS) + (R 2 (e), AG) 

- 7, {(£!, AS) + (e 2 , AG)} - ^ • {(£ 2 , VAS) - (£i, VAG)} 

- + bj {(e 2 ,A 2 S + 2a c AS)-( ei ,A 2 Q + 2a c AG)} - — ^ ^ F(e). 

It remains to extract the formally cubic term in e in the RHS (I3.14|) . We let: 



(3.14) 
(3.15) 



Gi(£) = Ri(e 



'' 1 \Q b \ p - 5 (pJ: 3 + 3SG 2 )e 2 - ^-i|Q b |P- 5 (S 3 + (p - 2)SG 2 )^ 



2 vr 2 

(p-i)|Q 6 r 5 (G 3 + (p-2)s 2 G) £l£2 , 



(3.16) 



G 2 (V) 



i? 2 (£) 



p — 1 



IQ b r 5 (e 3 + (p-2)Gs 2 )e 2 



p — 1 



|Q fe r 5 (pG 3 + 3GS 2 )e 2 



(P - l)|Q b | p - 5 (S 3 + (p - 2)Q 2 S)e l£2 , 



(3.17) 
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and eventually arrive at the following algebraic virial identity: 

"° ((f ' AS ) " (f - Ae ) " + + «-- A2 > " <- Ae »> - 

<%)<7 C + a c (2£(Q 6 ) + \Q b \\ 2 + t s |Q 6 || 2 ) - 2X 2( ~ 1 ~^E 
- {ex, Re{M> b )) - (e 2j Iro(A¥ 6 )) + flp(e, e) - 7s {(ei, AE) + (e 2 , AG)} (3.18) 
• {(e 2 , VAS) - (ei, VA9)} + E(e,e) + (Gi(e),AE) + (G 2 (e), AG) 

^ + &) {(e 2 ,A 2 S + 2 f T c AS)-(e 1 ,A 2 e + 2 CTc AG)} - | F(e), 

where the virial quadratic form H p can be expressed in the form: 

H p (e, e) = f |Ve| 2 + | „ • VQ p Q^ 2 e 2 



A 



(3.19) 



and the error E{e,e) simply comes from the error made by replacing Q b by Q in 
the potential terms and is easily estimated thanks to Proposition 12.61 bv: 

\E(e,e)\ < 5{p) (J |Ve| 2 + j \e\ 2 e-^ (3.20) 

with 5(p) — ► as p — > p c . 

step 2 Estimates of nonlinear terms and coercivity of the quadratic form. 

Observe from the construction of Q b that 



ae , _\ fas ._\ _ dQ b 

,AE — — ,AG = —Im AQ b , 



8b' J \db' J \ ' <% 

= -Im (AP b - ib^n, ^ - i^P^j = + 0(\b\ + <r c ).(3.21) 

We then inject the orthogonality conditions (lTMD . (l2~55l . (l2~56ll . (12371) and (l3~2TT) 
into (|3.18p to derive the algebraic identity: 

^^(l + <5(p))6 s (3.22) 
= 5{p)a c + a c {2E(Q b ) + |Q b | 2 2 + %\Q b \ 2 L2 ) - 2X 2 ^^E - (e 1} lte(Atf 6 )) - (e 2 , /m(A* 6 )) 
+ flp(e, E ) - 7s {(ei, AE) + (e 2 , AG)} - y • {(e 2 , VAS) - VAG)} 



+ E(e,e) + (Gi(e), AS) + (G 2 (e),A0) - [ F(e) 



with 5(p) — > as p — > p c . Let us now estimate all the terms in the RHS of (|3.22l) . 
First observe from (|2.34|) . (|2.36|) and the smallness of e the non degeneracy: 



a c \Q. 



|2 

S(p)a c + a c (2E(Q b ) + \Q b \ 2 L2 + y,\Q b \j?) > (3.23) 

The nonlinear terms are estimated as in [13], [15] using homogeneity estimates, 
Sobolev embeddings and the bootstrap bounds (|2.61l) . (I2.63P as for the proof of 
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(POTT) : 



E(e,e) + (Gi(e), AS) + (G 2 (e), AG) 



< <5(p) 



y | Ve i 2 + y | £ | 2 e -i«i) +r 



p+ i 



F(e) 



l+2 

6 



We now focus on the quadratic terms in (|3.22l) . We inject the estimates of the 
geometrical parameters (|3.7p . (13. 8j) . Using the simple bound 



e c|y| [|Q P -Q Pc | + |VQ p - VQ Pc 



as p — » p c 



we may view the obtained quadratic form together with H p as a perturbation of the 
one obtained in the L 2 critical case: 



H p (e,e)-%{(e 1 ,AZ) + (e 2 ,AQ)} 



A 



•{(e 2 ,VAE)-(e 1 ,VAe)} 



with 
H{e,e) 



" V£ " + N ( 1+ N 



Q P N C y • vQpA + jj I Q& v VQ P A 



-(euL+DiQpJie^DQpJ 



\DQ Pc \ L 2 

|^( £)£ )| <S(p) ^ |Ve| 2 + y | £ | 2 e-l^+r 



2-30^0 
b 



We now recall from |14j the following coercivity property which is a consequence of 
the spectral property stated in the introduction: Ve = £\ + ie 2 E H 1 , 



H[e,e) > co / |Ve|" + 



|e| 2 e-l"l 



- {(ei, Qp e ) 2 + (ei, \y?Q Pc ? + (ei,yQp e ) 2 + (e 2 , DQ Pc ) 2 + (e 2 ,D 2 Q Pc ) 2 + (e 2 ,VQ Pc ) 2 } 

for some universal constant Co > 0, and hence our choice of orthogonality conditions, 
(|1.22|) and the degeneracy estimates (|3.5I) . (|3.6p ensure: 



H(e,e)> C j[ l\Ve\ 2 + 



|el 2 e-W 



3 

r? 



(3.24) 



for Uq = vq(jp) > small enough. We now inject (|3.21l) . the orthogonality condition 
(|2.57j) . the nondegeneracy estimate (|3.23l) . the estimates on nonlinear terms (I3,24p 
and the coercivity property (I3,24p into (I3.22p to derive 



> 



co 



a c + 



\e\ 2 e-M 



1 b 



\Ve\ A + 

- |(ei,Jfe(Atf 6 )) + (£ 2 ,Jm(Att 6 ))|. 
A crude bound for the remaining linear term is derived from (|2.26|) (|2.30p fj2 . 3 1 [) : 



\(e 1 ,Re(A* b ) + (e 2 Jm(Ay b )\<r 1 - Cll + ^(a c + J |Ve| 2 + J \e\ 2 e^A, (3 



25) 



and (|3.12j) follows for r] < rj(vo) chosen small enough. 
This concludes the proof of Proposition [ 
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3.3. Refined local virial identity and introduction of the radiation. We 

now proceed through a refinement of the local virial estimate (13 . 12|) and adapt the 
analysis in |15j . 

We start with introducing a localized version of the radiation Cb introduced in 
Lemma 12.31 due to the non-L 2 slowly decaying tail (|2.22|) . Let a radial cut off 
function XA{f) = X (3) with x( r ) = 1 for < r < 1 and x(r) = for r > 2 with 
the choice: 



A = A(t) = e 2a ^) so that T b 2 < A < T b 2 , 

for some parameter a > small enough to be chosen later and which depends on r\. 
Let 

Cb = XA(b = Cl + <2- 
Cb still satisfies size estimates of Lemma |2~3"1 and is moreover in H 1 with estimate: 

2 



6(2) 



(3.26) 



l(i + |y|) 10 (IO.I + |vo ) |)|£ 2 + 



\10 



Ob 



+ 



V 



(Kb 

db 



< rl' Cv . (3.27) 



L 2 



From (I2.18p . the equation satisfied by Cfe is now: 



AO, - Cb + i^Cfe = *i 0) + i 7 



with 



(3.28) 



F = {A XA )Cb + 2V X A ■ VC b + • V XA Cb- 
We then consider the new profile and dispersion: 

Qb = Qb + Cb, e = s — Cfe 
and claim the following refined local virial estimate for e: 

Lemma 3.5 (Refined virial estimate for e). There holds for some universal constant 
c 2 > 0: 

{fi(s)} s >c 2 ' 



(3.29) 



with 



Ms) = --ImU yVQ b Q b ) - Im ( Cb, AQf> ) + (e 2 ,ACi) - (ei,AC 2 ) 
Proof of Lemma 13.51 



(3.30) 



The proof is parallel to the one of the local virial estimate (|3.12|) up to the es- 
timate of the remaining leading order liner term (|3.25|) for which will use a sharp 
flux computation based on (|2.22p . 



step 1 Estimate on ^>b 
Let 



^ b = -io- c Hb 
with explicitly: 



(2) 



dQb 
db 

dCb 



AQ b + Q b - ibAQ b - Q b \Q; 



iP-i 



-^o"c^ b — + ibo c Cb + * 
ob 



(i) 



(Q b + £b)\Qb + tb\ 



p-i 



i,( 2 ) 



(3.31) 



Qb\Qb 



.p-i 
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where ^ is given by (|2~2"9"1) . Observe from (l2~50j) . (EOT]) . (l3~2"T|) and the degeneracy 
on compact sets (|2.24l) that: 



\10 



\( 1 + \V\J U*6 
for some universal constants c, C > 



(l*i 2) l + |v^ 2) | 



Ii 2 <r fc 1+C + C7a c 2 



(3.32) 



step 2 Rerunning the local virial estimate bound. 

We now turn to the proof of (I3.29P and propose a small short cut with respect to 
the analysis in [15]. Let us indeed rerun exactly step 1 of the proof of Proposition 
3.31 with the new profile Q b an d variable e. Recall indeed that the whole algebra is 
completely intrinsic to the equation of the profile and relies only on the definition 
of respectively (|2.29|) . (I3.3ip . Moreover, the algebraic formula (|3.18|) does not rely 
on the choice of orthogonality conditions for e -which no longer hold for e— , and 
hence: 



~db 



~db 



,AQ 



£2, A 



~db 



+ ei,A 



oe 

~db 



+ 



|(£2, AE) - (ei,A6)| ( 



5(p)a c + a c (2E(Q b ) + |Q 6 | 2 2 + %\Q b \jp) - 2X 2 ^-^E 
(ea,i?e(A# 6 )) - (f 2 ,7m(A* 6 )) + H p (£,e) - % {(fi,AS) + (£ 2 ,AQ)} 

^ • {(£ 2 , VAS) - (fi, VAG)} + E(e,i) + AS) + (G 2 (e), AG) 



+ 6^) {(f 2 ,A 2 S + 2(T c AS) 



(fi,A 2 e + 2fj c AG)} 



p+ 1 



We first observe after an integration by parts that 
/ dO 



db ' 



AS 




Trim 

2db 



y ■ VQ b Q b -a, 



dO 
~db'' 




Using the orthogonality relation (I2.57p . we compute: 

yVQbSb) +(e 2 ,AS)-(£ 1 ,AQ) 



—Im 
2 



1 



Im( y- VQ b Q b )-2Im( ( b , AQ b + (e 2 , ACi) - (ei, A£ 



= /i- 

Next, all the terms in (|3.33l) are treated like for the proof of (13.12p except the linear 
term involving ^ b . Arguing as in [15] and using in particular the degeneracy of 
£b on compact sets (|2.24p to treat the scalar products terms in e, we arrive at the 
following preliminary estimate: 



{/l} s > c 



|W| 2 + 



\i\ 2 e 



\v\ 



CX 2{ - l ' ac) E n - Ca r 



1 b 



- (£i J i?e(A* 6 ))-(£ 2 ,Jm(A* 6 )) (3.34) 

for some universal constants C, zq > and with /i given by (|3.30p . It remains to 
estimate the leading order linear term. We first estimate from (I3.3ip . (|3.32p : 

-(£ 1 ,i?e(A^ b ))-(e 2) /m(A# 6 )) > (i u Re(DF) + (e 2 , Im(DF)) (3.35) 



10 



iVel 2 + 



\i\ 2 e-M 



.l + CQ 



(3.33) 
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To estimate the remaining linear term, we proceed as in [15] and split i = e — (f. 
(£i,#e(I}F))+(e 2 ,Im(£F)) = ^^ 

The last term is the flux term for which the following lower bound can be derived 
from (l2~22D : 

- (Ci, Re(DF)) - (C 2j Im(DF)) > c T b . (3.36) 
The other term is estimated from Cauchy Schwarz and a sharp estimate on F from 
M and (|222) , (1223): 



\(e l ,Re(DF)) + (e 2 ,Im(DF))\ < (j kN (/ l*f) 

We refer to [TS], step 4 of the proof of Lemma 6, for a detailed proof of (|3.36|) . 
(POT!) . Injecting (EOT)]) . ([3~3"Tjl into (l3~3"5D and (l3~3"4l now yields (l3~2HD . 
This concludes the proof of Lemma [ 



3.4. Computation of the I? flux. We now turn to the computation of 1? fluxes 
which are the key to get upper bounds on the far away localized L 2 term which 
appears in the RHS of (I3,29p . The obtained identity displays new features with 
respect to the analysis in [15] which reflect the I? super critical nature of the 
problem. We introduce a radial non negative cut off function 4>{r) such that <p(r) = 
for r < i, cf)(r) = 1 for r > 3, \ < <p'(r) < \ pour 1 < r < 2, cf)'(r) > 0. We then let 

<l>A(s,r) = (f> 

with A(s) given by ([3726]) . 



A(s) 



Lemma 3.6 (L 2 fluxes). There holds for some universal constant Cs,zq > and 
s > 0: 



1 



r2A 

A 2CTc / </> A \e\ 2 \ >c 3 b I \e\ 2 -Tl +Z0 -Tg I \Ve\ 2 . (3.38 



\2* c 

Proof of Lemma 13.61 



.s 



A 



Take a smooth cut off function x(i, x). We integrate by parts on (II. 1|) to compute 
the flux of 1? norm: 

x(t, x)\u(t, x)\ 2 dx 1 = — [ dtx{t,x)\u(t,x)\ 2 dx + Im (Vx • Vuu) . 



2 1./ I * 2 . 

We apply this with x(^ x ) = &a( ^t$- ) and inject the decomposition (|2.53p . Recall 

by construction that Qb(y) = c c 7fe(y) for |y| > | which is uniformly exponentially 
decreasing and hence using (I2,58p . its contribution near A generates terms which 
are negligible with respect to the leading order Tfe. We get after a bit of algebra 
using also the bootstrap estimates: 



s : + b)Jv v^kl 2 - y± ■ J vMe? - rl +zo (3.39) 
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for some universal constant zq > 0. Prom the choice of 

1 



10 



\e\ 2 > 



e\ 2 > 



10 



/ 1 y_ 

A 



\e\ 2 > 



40 



2A 



\e\ , (3.40) 



and also from the choice of A and the control of the geometrical parameters: 



Moreover: 



y ■ V0A|e| 



+ 



ds 



A, ,2 



~ 1000 



(3.41) 



< 



< 



1 

I 

40 
6l2 



|Vd 



l VE l + 40 



40 



< ^ ^ 4 kl 2 + r 



.4 



a) |£| 



|Ve| 



(3.42) 



Injecting (I3,40p . (|3.41l) and (|3.42l) into (|3.39j) yields (|3.38j) and concludes the proof 
of Lemma 13. ( 



3.5. L 2 conservation law and second monotonicity formula. We now couple 
the estimates (I3.29p and (|3.38p together with the L 2 conservation law to derive 
a new monotonicity formula which completes the dynamical information given by 
(pU2l) . 

Proposition 3.7 (Second monotonicity formula). There holds for some universal 
constant C4 > 0: 



{J} s > c 4 b (r b + j |Ve| 2 + J \e\ 2 e-\y\ 



-Or. 



C 4 



(3.43) 



with 



J{s) = (J \Q b \ 2 - Jq 2 ) +2(ei,E) + 2(e aj 0) + J \l- <j> A )\e\ 2 (3.44) 
- c 3 c 2 (bh(b) - f h{v)dv + b{(e 2 , ACi) - (e^ACa)}! , 



where 03,02 are the universal small constants involved in fl3.29[) . (|3.38l) . and: 

1 



h(b) = -Im{ j yVQ b Q b ) + (0, AS) - (£, AG). 



(3.45) 



Proof of Proposition 13.71 

step 1 Coupling 13331) and (l3~38l) . 

Let us multiply (|3.29p by bci'- 



be 2 



J |Vef + j lelVM+lO < b{c 2 fi} s + b £ A \e\ 2 + bo, 



(3.46) 
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We then integrate by parts in time using (|3.30l) : 

b{c 2 fi} s = cfc{&/i(&)- jf f 1 (v)dv + b(E 2 ,ACi)-b(s 1 ,AC 2 ) 

- c 2 6 s {(e 2 ,ACi) - (ei.ACa)} 
and estimate from (|3.7|) and fl3.2Tj) : 

c 2 6 s {( £2 ,ACi)-(£i,AC2)} 

c 2 6 s {(£ 2 ,ACi) - (ei,AC 2 )}| + |c 2 6 s {(Ci,AC 2 ) - (C 2 ,ACi)} 

rl +Z0 + b -§(J l^l 2 + / l^i + r fc ). 

Injecting this into (|3.46|) yields: 

2 'I |W|2 + / |f|2e "' yl+rfe ) " C2 { 

+ b 



< 
< 



be? 



6/1(6) 
2.4 



(3.47) 



/i(t0d« + &(e3,ACi)-&(e 1 ,A<2) 



\e\ + ba c . 



We now inject the control of 1? fluxes (13.38P and obtain for a > Crj: 



|Vef+ / |f| 2 e -|yl +r 6 ) < c 3 c 2 {bf 1 (b) 



+ i{ A2CTC / ^N 2 } s + 6c 3( r c . 
step 2 Injection of the L 2 conservation law. 
We now rewrite the L 2 conservation law as follows: 



/i(v)cfo + &(e2 ) A6 ! )-6(ei,ACi) 

(3.48 



K| 2 = A 



2<r r 



\Q b \ z + 2Re(e,Q b )+ / |e 



which yields: 



A 2<x c 



i2cr c 



1 2(T C 



4>a\e\ 2 

(l-</> A )M 2 + 2i?e(e,Q^) + ^ \Q b 



(l-cf> A )\£\ 2 +2Re(e,Q b )+ / |Q. 



(l-^)k| 2 + 2 J Re(e,Q fc )+ / \Q 



Now the Hardy type bound 

j(l-d> A )\s\ 2 <CA 3 (| |V,| 2 + J \e\ 2 e-\y^j 

together with the choice of A (|3.26p . the bootstrap bound (|2.6ip and the control 
(|3.7p yield the rough bound: 

A, " 



(3.49) 



Or 



A 



(l-0 A )kr + 2Re(e,Q b )+ / \Q 



< Cboc. 
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We may thus rewrite (I3.48P as: 

2 (J | W | 2 + J | £ f e -M +r ^ <Cba c 



bc 3 c% 



+ I- (J \Qb\ 2 ~ jQl\-2Re(s,Q b )- |(1-<^)| £ | 2 



+ c 3 c 2 



b 



6/1(6) " / h(v)dv + 6(e 2 , AC2) - 6(ei, ACi] 







which is (|3.43|) . This concludes the proof of Proposition l3~77l 

4. Existence and stability of the self similar regime 

This section is devoted to the proof of the main Theorem 11.21 We first show 
how the coupling of the monotonicity formulae (|3.12p . (I3.43P implies a dynamical 
trapping of b and a uniform bound on e which allows us to close the bootstrap Propo- 
sition 12.121 We then conclude the proof of Theorem 11.21 as a simple consequence of 
these uniform bounds. 

4.1. Closing the bootstrap. We are now in position to close the bootstrap and 
conclude the proof of Proposition 12.121 

Proof of Proposition 12.121 

step 1 Pointwise bound on e. 

Let us start with the proof of the pointwise bound on e (|2.67|) . We argue by 
contradiction and assume that there exists s 2 G [so>s*] such that: 

J \Ve(s 2 )\ 2 + J \e(s 2 )\ 2 e-^>Tl~^. 

A simple continuity argument based on the initialization of the bootstrap estimate 
(|2.52|) implies that there exists [53,54] C [so,s*] such that: 

J |Ve( S3 )| 2 + J \e(s 3 )fe-M = 1^°, J \Ve( Si )\ 2 + J | £ ( S4 )|VW = I*"** , 

(4.1) 

and 

V, G [s 3 ,s,], J \Ve(s)\ 2 + J \s(s)\ 2 e~\y\ > T^. (4.2) 
From (|4.2|) and the first virial monotonicity ()3.12|) . we have: Vs G [s 3 , S4], 

6 s >c 1 (r 6 1 - 7 ^-r 6 1 - l/ ° 2 )>o 

and hence 

6(s 4 ) > 6(s 3 ). (4.3) 
On the other hand, using the lower bound 

j | Ve f + j \i\ 2 e-\y\ >UJ \Ve\ 2 + j \e\ 2 e^A - T^' 

together with (|4.2|) . (|2.58|) and the second monotonicity formula (|3.43p . there holds: 
Vs G [33,34], 
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and hence 

J{sa) < J (s 3 ). (4.4) 
We now claim the following upper and lower control of J: 

J{S) h(b(s)) { < CA 3 (J | V£ |2 + J | £ |2 e -|„|) + T l-Ca + C(Jc) (4-5) 

where /2 given by 

f 2 (b,a c ) = (J \Q b \ 2 - J (A - c 3 c 2 (bh(b) - J fi(v)dv} (4.6) 

satisfies 

V6*>6 2 >6i, / 2 (6 2 ) > / 2 (6i) - C(7 C , ^-Ca c </2(6i)<C(&? + ff c ). (4.7) 
Let us assume (jOj) . ([4~7D . Then (j4~4) . ([43]) imply: 

h(b(s,)) - + I (/ |Ve( S4 )| 2 + | \s( Si )\ 2 e-^ - Ca c < J( S4 ) < J(s 3 ) 

< f2(b(s 3 ))+CA 3 (s 3 ) (| \Ve(s 3 )\ 2 + J | £ ( S3 )| 2 e-^l) + T\$* + Cv e 
and hence from the monotonicity (|4.3I) . (j4.7|) and the controls (|4.1|) . (j2.58|) : 



7uo-Ca + c < pl-S^o 

6(53) — £>(S4j 



for a = C?7 > and 77 chosen small enough, a contradiction which concludes the 
proof of (1227}. 

Proof of (|4.5p . (|4.7p : It is a standard consequence of the coercivity of the linearized 
energy with our choice of orthogonality conditions, |15| . Indeed, we rewrite J given 
by (|3.44l) using the conservation of energy (13.9(1 and the orthogonality condition 
(12371) : 

J = f 2 (b,a c )+2\ei,Re(^ b )-a c fx b ^\ + 2 (e 2 , Im(V b ) + a c fi b ^ 

- c 3 c 2 b{(s 2 , ACi) - (ei, AC2)} + 2E(Q b ) - 2X 2 ^-^E 
+ (M + (e), £l ) + (M_(e),e 2 ) - | </> A k| 2 - ^ | F(e). 

The upper bound in (|4.5p now follows from the Hardy bound (|3.49l) and the degen- 
eracy (|2.36|) . For the lower bound, we recall the following coercivity of the linearized 
energy which holds true for A large enough i.e. \b\ < b* small enough: 

(M+(e),ei) + (M-(e),e 2 )- J 4> A \e\ 2 > c 3 (J |Ve| 2 + J \s\ 2 e^A 

- {(ei,Q Pc ) 2 + (ei, \y\ 2 Q Pc ? + (ei,yQ Pc ) 2 + (e 2 ,D 2 Q Pc ) 2 } , 
^3 



for some universal constant c 3 > 0, see Appendix D in [15]. The choice of orthog- 
onality conditions together with the degeneracy (|3.5p now yield (|4.5p . (|4.7|) is now 
a direct consequence of (|2.34|) . (|2.35p . Indeed, first observe from (|3.45|) . (|2.33p and 
the smallness of the radiation given by Lemma l2~3l that : 

dMxib) 



/ x (6) = Mi (b) + 0(a c ) with Mi(0) = and 



db 



C{p) 



31 

for some universal constant C(p) > 0. Hence (|2.34ft . (|2.35ft imply: 

f2(b,a c ) = (J \Q b \ 2 -jQl]- c 3 c 2 (bh(b) - J h(v)dv 
= M{b) + 0{a c ) 



with 



dM(b) , . , . „ , b , , 
7 > b(c (p)-Cc 3 c 2 ) > -co{p c ) 



db ~ " 2 

provided the constants C2, C3 in (I3.29p . (I3.38P have been chosen small enough, and 
the monotonicity (|4.7ft follows. 

step 2 Dynamical trapping of b. 

We now turn to the core of the argument which is the dynamical trapping of b 
(222}. We recall from (l2~22l that 

Assume that there exists S5 G [so,s*] such that a c > T b ,?, then from (|2.49|) and 
a simple continuity argument, consider sq G [so^s) the largest time such that 

/ g 0(2) \ l ~ U 

a c = I e b ( s e) J , then b s (sQ) < by construction while from (13. 12ft : 

6 s ( S6 )> Cl (a c -r^)> Cl ^e- 2 ^y V °-Tl( s fj >0 
and a contradiction follows. 

Similarly, assume that there exists sj G [so,s*] such that 

From (|3.44ft . the pointwise bounds (I2,58p . (12.611) . and the value of the L 2 norm of 
Q b given by J23H), l[23H11 . there holds: 

l^-^ 2 |<^ 2 (4-9) 
for some universal constant do > 0. Hence (|2.49p . (14. 8ft imply: 

rfo / V V d o 



We then consider -using again (|2.22p - the largest time ss G [sol's?] such that a c = 

2 9(2) \ i+^O 

e I ( then (J%(s 8 ) > by definition while from j3~43]), (US), J2Z22]): 



{ J} s (as) > 6(ss) 



>0, 



and a contradiction follows. This concludes the proof of (|2.64ft . 
step 3 Control of the scaling parameter. 
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We now turn to the control of the scaling parameter X(t). From (|2.49l) . the 
dynamical trapping of b (I2,64p implies: 

Vt€[0,T*), (1 - ^)b < b(t) < (1 + v%)b . (4.10) 

Hence the upper bound (|3.7p . the control (|2.67l) and (|4.10p ensure: 

Vi€[0,T*), 0<(l-2^ 3 )6 <-y = -A t A<(l + 2t, 3 )6 . (4.11) 

In particular, A is nonincreasing while b(t) is trapped from (I4.10P and thus (|2.50|) . 
(12311 imply (I2~65ll . (P2~66l) . 

step 4 Control of the solution in H a . 

It remains to close the H a estimate (|2.68p which is the key to the control of the 
nonlinear term (|2.63p . We use here the fact that the blow up is self similar and 
strictly H 1 subcritical so that 

N N N 2 



" 2 p+1 2 p-1 

In other words, norms above scaling can be controlled dynamically in the bootstrap 
as was for example observed in [21], [23J, |22j. 

It is more convenient here to work in original variables. Consider the decomposition 

u(t, x) = Q sing (t, x) + u{t, x) = (Q b + s) (t, X ~, X ^ ) e^ l) , 

A^(i) V Mi) ) 

then first observe by rescaling and the trapping of b (|4.10l) that (|2.68|) is implied by: 

pi— 45^o 

To prove (|4.12p . we write down the equation for u and use standard Strichartz 
estimates, see |T], for the linear Schrodinger flow. Indeed, the equation for u is: 



idtu + Am = —6 — f(u) 



with: 



£ — idtQ sing ~\~ AQsing "f" Q sing\Q sing\^ 
1 



1 

A 2+ ^ 



ib s ^ + AQ b -Q b + Q b \Q b \ p - 1 - i^AQ b - & ■ VQ b - %Q b 
ab A A 

dQb f X s \ x s 
+ i(b s - a c ^ h )-^r "My + H ^Qb - «y • VQb - IsQb 



x 

t, 



and 



/('^) — {Qsing ~i~ ^)\Q sing ~i~ Qsing | Qsing |^ • (4'13) 

Let £ G [0, T*), we write down the Duhamel formula on [0, t]. Following [T], we 
consider the Strichartz pair: 



i\T + a(p-l)' (p-l)(JV-2<r)' 7 2 

and estimate from Strichartz estimates: 

nvr*ui L oo ^ < iivrcioi^ + iivrari t2 + iivr/(«)i r v r ,. (4.15) 
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Using the fact that A is nonincreasing together with (12, 52ft . we obtain: 

,s(1-md) 



^^^vk^- (416) 

We claim: 

r £(i-i5no) 

w £ h.^. £ m^>' (4 - 17) 

r ±(i-4i^ ) 

IIVr/W«))lL» S plp^ («8) 
which together with (jilUjl . ([47T5]) . ([4TBD yield P~T2j) . 

Proof of (|4.17|) : From the estimates on the geometrical parameters (13, 7|) . (|3.8p . the 
degeneracy estimates (|2.30p . (12.3ip . the pointwise bounds (12.67P and (I2,64p . and the 
freezing of b (14TTO . there holds: Vi 6 [0,t*] s 



ivrf(t)i L2 < 



i 



[A(t)]2+(«r-«Tc) 



-* 6 + i(b s - o-cUb)—- - i I + b ) AQ b - i-f- ■ VQ b - ^ s Q b 



db V A / A 



i 

J / / ,0 / , , _U,I -1-11;, \ 2 



[A(t)]2+(<r- CTc ) 

r |(i-iawo) 



< 7TTTTTT77Z — 7T ( / I V £ | 2 + / + T^" 1 



- [A(t)] 2 +(--^) ' (4 ' 19) 
We now observe from the self similar blow up speed estimate (|4.11l) : Vq > 2, 
dr C /■* A t C . 



[A(r)]« " 6o7o [A(r)]*-i " (g - 2)6 [A(t)]^ " 
Integrating (I4.19P in time and using (|4.20l) yields (14,l?p . 

Proof of (I4.18P : This estimate follows m £/ie bootstrap using the fact that the blow 
up is self similar and that u is small in H a for a > o c after renormalization. Indeed, 
we first claim from standard nonlinear estimates in Besov spaces: 

\M°f(u)\ Lr ,<^^\\Vfs\ L *, (4.21) 

where a is defined by: 

N N 2 , 

a = a+ — =a + -. 4.22 

2 r 7 

The proof of the estimate (|4.21|) is postponed to the appendix. From direct check, 
a < a < 1 providing a has been chosen close enough to the critical scaling exponent 
a c himself close enough to 0. We may thus interpolate between a and 1 and use 
(12T621) . (12T671) and (ET221) to estimate: 

\M a e\ L 2<\\V\°e\lr\Ve\lr <T^ V J _ ( 4 . 23) 

Provided p is chosen close enough to p c , and a is chosen close enough to 0, we obtain 
from and (ICTll . (l4~23D : 

HVl^k^rf- 40 ^. (4.24) 
Injecting this into (|4.21|) yields: 
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Now from (fl~T4) . (14T221) . there holds: 

2 2n / 1 \ v + 1 

(or - <r c )p - — = p(a - a c ) + - 2 1 - - = p{a - a c ) - 2 + 2 " 



/ T V T / T 

JV- 2a" 



(O- - <7 C ) + (p - 1) 



cr - <T C + 



= CJ - CJ C . (4.26) 

In particular, (a — cr c )pj' = 2 + 7' (a — cr c ) > 2, and we may thus inject (|4.20p into 
(|4.25l) to conclude: 

£(l-41wo) 

for z/o > small enough thanks to p > 1, this is (|4.18p . 

This concludes the proof of the bootstrap Proposition 12.121 
4.2. Proof of Theorem 11.21 We are now in position to prove Theorem 11.21 

Proof of Theorem [TT2l 

Pick vq > 0, p G (Pc,P*(vo)) and uq G O corresponding to 60 = b*(p) as given by 
Definition 12.81 Note that (|1.12h follows from (|2.49|) . Let u(t) be the corresponding 
solution to (jl.ip with maximum lifetime interval on the right [0, T), then from 
Proposition 12.121 u(t) admits on [0, T) a geometrical decomposition 

u(t,z) = -^—(Q m + e) (t, ^P-) e i7(t) 

which satisfies the estimates of Proposition 12.121 This implies in particular (|1.14|) . 

step 1 Finite time blow and self similar blow up speed. 

Recall ffiTD : 

V* G [0, T), (1 - v$)b < -X t X < (1 + Re- 
integrating this in time first from to t yields: 

1 A 2 

Vt€[0,T), (1 - ^ 2 )M < 77A 2 , and hence T< 

so that the solution blows up in finite time. From the H 1 Cauchy theory, \S/u(t)\ L 2 — > 
+00 as i — > T and hence from (|2.67j) . A(i) — > as t — > T. We thus integrate (14. lip 
from t to T to get: 

Vt G [0, T], (1 - ^ 2 )6o(T - i) < ^ < (1 + ^ 2 )6 (T - t) 
which implies (|1.16|) . 

step 2 Convergence of the blow up point. 

From (13. 8|) and Proposition 12.121 we have the rough bound: 



35 



and thus from (|1.16l) : 



1 

T r-T p4 



o Jo ^b {T-t) 



and (|1.15h follows. We moreover get the convergence rate 
z(t) - x(T) 



i 



X(t) 



f T dr 
~ J ^/2b (T-t) Jt ^2b (T - t 



< t= = „ / j= i <K- ( 4 - 28 ) 



step 3 Strong convergence in H s for < s < a c . 

We now turn to the proof of ()1.17p . Pick < s < a c . Let < r <C 1 and 

< t < T — r, let u T (t) = u{t + r) and v(t) = u T (t) — u(t), then v satisfies: 

iv t + Av = ulu^" 1 - it T |it T | p_1 . (4.29) 
Consider the Strichartz pair 

N(p+1) 4(p+l) 2 JV JV 



N + s(p-iy 7 (p-l)(iV-2s)' 7 2 r' 
then from standard nonlinear estimates in Sobolev spaces -[I]-, we have: 

||V| S [n^l^ 1 - n.ln.r 1 ] \ Lr , < \\V\ s u\ p Lr + ||V|V|^ < ||V| ff u|^ + ||V|V|£ 2 
with 

JV Af 2 

a = s + — = s + -. 

2 r 7 

Now observe that a — > jM^ > as p — > p c and hence a c < a < a < 1 from (|2.48|) for 
p close enough to p c . We thus estimate from the geometrical decomposition (|2.53p 
and the bounds (12371) . (l2~68l) : 

l|Vr-u| L 2 < _J_||V| ff (Q 6 + e)lxa < t^— • 
We conclude using (|4.26l) : 

nvr H -r' - «rKM i <t _^, s (f nvr<;) * < (f p^^) 

~ (f [A( T )A<~.> ) y ~°»'~ 1 * 4 - 



from the scaling law X(t) ~ y / 2bo(T — t) and s — <7 C < 0. By running the standard 
Strichartz estimates -[!]- on (|4.29p . we conclude that: 

( f' T dr \^ 

I|vN^ t _ t)L2 <||v|M0I, 2 + (/ [A(T)]2+y(a - ge) j » 

and the continuity ii 6 C([0, T), iT s ) now implies that u(t) is Cauchy sequence in 
H* as t -> T, and Q7I) follows. 

Remark 4.1. iVoie £/iai i/ie case s = a c in (|4.30l) /mrfs i/ie logarithmic upper 
bound on the critical norm (|1.20p . 
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step 4 Behavior of u* on the blow up point. 

It remains to prove fjl . 18[) which follows by adapting the argument in [16J. 
Let a smooth radially symmetric cut off function x( r ) = 1 f° r r < 1 and x{ r ) = 
for r > 2. Fix t £ [0, T) and let 



with Aq given by 



R(t) = A X(t) 



a 2 ^ 
Aq = e b o 



(4.31) 
(4.32) 



We then compute the flux of L 2 norm: 



d_ 

7h 



X 



x — x(T) 

R{t) 



R(ty 



Im[ Vx 



x — x(T) 



R(t) 

< —Hr)fU< ' ' 



Vu(t)u(t) 



R(ty v n Hi ~ R(t) [A(r)] 1 - 2 ^ 



where we used (I2~67ll . (I2~68ll . We integrate this from t to T, divide by R 2cTc (t) and 
get from (fTTT7l) : 



1 



X 



< 



< 



T 



x — x(T) 

W) 

d.T 



,*|2 



R 2 ^(t) 



X 



R 2a -{t) 



x — x{T) 

W) 

T 



\u(t)f 



< 



(It 



1 



< 



1 



K c+ %' A 



1/2+2<t c 



(4.33) 



where we used the self similar speed ()1.16|) and (|4.32p . On the other hand, we have 
from (12331 : 



R 2 °c(t) 



X 



x — x(T) 

R(t) 



|u(*)| s 



R 2a -(t) 
1 



A 



2o> 



X 



X 

I 

T 



X(t) ( x{t)-x{T) 

y + 



R(t) 

y + 



X(t) 
x(t) - x(T) 



\Q b + e\ 2 (y)dy 



\Q b + e\ 2 (y)dy. 



Now observe from the Hardy type bound (|3.49p . (|4.32|) and the bound (|2.67p that: 



/ 

J\v\ 



\y\<WA 
and hence fl4.28j) ensures: 



X 



x(t) - x{T) 

W) 



\Qb + e\ 



\Q p \ 2 (l + 5(p)) 



with 5(p) — ► as p — > p c . Injecting this into (|4.33p yields: 



R 2 ^(t) 



, x-x{T) \ 2 

v 1 ~1W 



A 



2(7 c 



\Q p \ 2 (l + 5(p)) + 



1 



1/2+2<t c 



A, 



o 



We now let t ^ T ie R(t) -> from (|4T3H and (fTTHll follows. 
This concludes the proof of Theorem 11.21 
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Appendix 

This appendix is dedicated to the proof of (14.2ip . The authors are grateful to 
Fabrice Planchon for having put them on the right tracks. 
After rescaling, we have: 

\M a f(u)\ Lr/ =_^_ y ||vr(F(Q 6 + e )-F(e))| Lr /, 
where the function F : C — > C is defined by: 



F(z) = \z\ v ~ x z. 



(4.34) 

(4.35) 
(4.36) 



Thus, (|4.21l) is equivalent to: 

\\Vr(F(Q b + e)-F(e))\ Lr , <\\Vfe\ L 2. 
We now concentrate on proving (|4.36p . We first rewrite F(Q b + e) — F(e) as: 

F(Q b + e) - F(e) = (J d z F(Q b + re)d^j e + (j d z F{Q b + T e)dr^j e. (4.37) 

Both terms in the right-hand side of (I4.37P are treated in the same way. Thus, for 
simplicity we may only consider the first term in the right-hand side of (|4.37j) . We 
introduce the real number q such that 

111 , , 

- = ---. 4.38 

q r r 

In particular, using the definition of r (I4.14p and usual Sobolev embeddings, we 
have for any function h: 

\h\ L i^ )q <\M a h\ L r. (4.39) 
Also, using again Sobolev embeddings together with the definition of a implies: 

\M a h\ L r<\\Vfh\ L 2. (4.40) 

We also introduce the real number u > 1 defined by: 

1 1 1 



v r' (p — l)q 

()4.38|) and (|4.4ip together with standard commutator estimates -see [8]- yield 



(4.41) 



|V| CT (sj d z F(Q b + rs)dT S ] 



< livrd 



+ I £ Il(p- 1 )8 


ivr 











/ O z F(Q b + TE)dT 

Jo 
1 

d z F(Q b + TE)dT 



L'l 



(4.42) 



v 



which together with (I4.39P and (|4.40l) yields: 



Vr e / d z F(Q b + re)dr 



<\M a e\vlj \d z F{Q b + re)\ Lq dr 
+ J o \\V\ a [d z F(Q b + Te)}\ Lv dr\ 



(4.43) 



Thus, in view of (I4.37|) and (|4.43p . proving (|4.36|) is equivalent to proving the 
following bound: 

f \d z F(Q b + rs)\ Lq dr + C ||Vr [d z F(Q b + tb)]\ lv dr < 1. (4.44) 
Jo Jo 

Now, we have by homogeneity: 

Vr G [0, 1], \d z F(Q b + ts)\< IQ^" 1 + \e\^ 1 
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which together with tl4T24|) . (f4739)) and (OPT) yields: 



/ \d z F{Q b + re)\ Lq dr <l {\Q b \^l_ x) + \e\^ n )dT 

<f(i + iW)* (4 - 45) 

Thus, we have reduced the proof of (|4.2ip to the proof of the following bound: 

f 1 ||Vr [d z F(Q b + re)]\ LV dr < 1, (4.46) 
J o 

where f is defined in (|4.41h . 

From now on, we concentrate on proving (I4.46p . To ease the notations, we define: 

h T = Q b + re, <r < 1. (4.47) 

Recall -see p]- the equivalent expression of homogeneous Besov norms: VO < a < 1, 
VI < g < +oo, 



f CT sup \u(- - y) - u(-)\li 

\y\<t 



i 

2 \ 2 

dt 



(4.48) 



Also, recall that ||V| CT «|i8 < |u| 

We start by proving (|4.46p in the case 1 < N < 4. In this case, p > 2 since 
p > p c . Using the homogeneity estimate: 

Vu,v, \d z F(u)-d z F(v)\<\u-v\(\u\P- 2 + \v\P- 2 ), 

and the relation from (l4T38l) - ([4T4lD : 

11 1 1 p-2 

+ 



v r' (p—l)q r q{p — 1) : 
we first estimate from Holder and (|4.39p . (|4.40p : 

\d z F{h T ){--y)-d z F(h T )(-)\ LV 
\{h T {- -y)- h T (.))[\h T (- - y )r 2 + \K(-)r 2 ]\ Ll/ 



< 

< \h T (--y)-h T (-)\ Lr \h T \ p ~^ 1) 

< \h T (--y)-h T (-)\ Lr \\V\ a h T f- 2 



and hence from (|4.40D and (I4.48P : 

iivr \d z F{h T )\ \u < Hvr^uHivr^i^ 2 < iwfKf- 1 . (4.49) 

(I4T24D and the definition of h T (l4~4Tll yield: 

\\Vfh T \ L 2 < \\VfQ b \ L2 + \\Vfe\ L2 < 1, (4.50) 

which together with (|4.49p implies the wanted estimate (|4.46l) . 

We turn to the proof of (|4.46|) in the remaining case N = 5. In this case, p > 9/5. 
We define the real number 8 by: 

1 1 p-9/5 



6 v q(p — 1) 
Using the homogeneity estimate: 

V«,u, \d z F{u) - d z F{v)\ < \u - v\l (\u\ p -% + \v\ p -%), 



(4.51) 
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we first estimate from Holder, ([4391) . (f440l) and (1430]) : 

\d z F(h T )(--y)-d z F(h T )(-)\ L , 



< 




-y)- 


< 


M-- 


-y)- 


< 




-y)- 


< 




-y)- 



4 

49 

K T 
'If 



(4.52) 



We decompose the integral in fl4.48j> in t > 1 and i < 1. For t > 1, we use: 



M--v)-M-)l%fl <i + M--y)-MOI r ¥ 



so that: 




t CT SUp |/l T (- - y) - MOLiii 



IT 




(4.53) 



By usual Sobolev embeddings, we have: 

\hr\ B « < 



(4.54) 



where a\ is defined by 

1 o\ _ 5 a 
2~ ~5 ~ 49 ~ 5' 
Using the definition of r, 7, q, u, 9 and a±, we obtain: 

2 3d 

CTl = ~ + — 

7 4 

which satisfies cr < <7i < 1 for p close enough to p c and a small enough. Thus, the 
bootstrap assumptions (|2.61l) (I2.62p and the definition of h (j4.47j) yield: 



iivr i / lr | L2 <nvr i Q b | L2 + iivr i e| L2 <i 

which together with (|4.53p and (|4.54l) implies: 

2 




t-° su V \h T {--y)-h T {-)\\ e 

\y\<t 

For t < 1, we use: 

M- - y) - MOI^e < ^ + \hr{- - y) - MOI 

so that: 



(4.55) 



(4.56) 



LIT 
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By usual Sobolev embeddings, we have: 



\K\ .flf < \\VpK\ L 2, 



(4.58) 



'40 



where 02 is defined by 

1 02 _ _5_ 9a; 

2 ~ T ~ 49 ~ 20' 
Using the definition of r, 7, q, v, 9 and 02, we obtain: 

2 

0-2 = - + la 

7 

which satisfies a < 01 < 1 for p close enough to p c and <r small enough. Thus, the 
bootstrap assumptions (|2.61|) (I2,62p and the definition of h (I4.47P yield: 

l|Vr 2 /i r | L2 < \\VpQ b \ L s + ||Vr 2 e| L2 < 1, (4.59) 

which together with (I4.57P and (|4.58j) implies: 

2 



t a sup \h T {- -y) -h T (-)\\ e 

\y\<t LIT 



*| I. 

t 



(4.60) 



Finally, (OHll . dOSll . TOH) and (001 yield: 

HVn^F^)]^ <1, (4.61) 
which implies the wanted estimate (|4.46[) . This concludes the proof of (|4.2ip . 
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